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CALCULATIONS OF THE FLOW OF NATURAL GAS 
THROUGH CRITICAL FLOW NOZZLES 

by Robert C. Johnson 
Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio 44135 

ABSTRACT 

The mass flow rafce of methane and nineteen natural gas mixtures through 
critical flow nozzles has been calculated. The calculation assumes the flow 
to be one -dimensional and isentropic. The pressure range is 0 to 1000 pounds 
per square inch and the temperature range is from 450 to 700 degrees Rankine. 

From a study of the results, a simple empirical method for making this 
mass flow rate calculation is proposed. This method would apply to natural 
gas mixtures whose composition is known and whose components have no more 
than four carbon atoms. 

INTRODUCTION 

When critical-flow nozzles are used for metering the mass flow rate of 
natural gas, the conventional isentropic flow relations do not apply. These 
equations only apply to a perfect gas. A perfect gas is defined as one 
having a compressibility factor of unity and an invariant specific heat. 

Gases such as air and nitrogen closely approximate this ideal condition at 
room temperatures and pressures up to a few atmospheres. Natural gas, on 
the other hand, cannot be considered perfect even at pressures less than 
atmospheric because natural gas has an appreciable specific -heat variation 
with temperature. At higher pressures, the compressibility-factor variation also 
becomes important. 
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A number cf methods have been used to calculate the isentropic mass 
flow rate of nonperfect gases. One of these methods is that used in 
reference 1. In this reference, an isentropic exponent is calculated, 
using the thermodynamic properties of a real gas. This exponent is then 
used in the conventional critical-flow equation to calculate the mass 
flow rate. Reference 2, using the tabulated data of reference 3, presents 
tables of compressible flow functions for the one-dimensional flow of air. 
Reference 4, using the state equations of reference 3, presents tables of a 
critical flow factor for air, nitrogen, oxygen, normal hydrogen, parahydrogen, 
and steam that permit the calculation of the isentropic mass flow rate of 
these gases through critical-flow nozzles. Reference 4 refers to a number 
of other reports that describe other methods for making this type of calculation. 
Reference 5 presents the critical flow factor for nitrogen and helium, and, 
in addition, describes a computer routine for making these real-gas 
computations . 

In this report, the critical flow factors for methane and two typical 
natural-gas mixtures are presented. In addition, from a study of nineteen 
natural gases, an empirical method for calculating both the mass flow rate 
and upstream volume flow rate of natural gas through critical-flow nozzles 
is given. Biis method requires knowledge of the upstream pressure, temperature, 
and composition of the natural gas. These calculations encompass a temperature 
range of 450°R to 700 °R and a pressure range of 0 to 1000 psia. 

NOMENCLATURE 

A q Constant in state equation (equation 13 ) 

A^ Area of nozzle throat 

a Constant in state equation (equation 13 ) 
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Constant in equation 23 

Constant in equation 24 

Constant in state equation (equation 13) 

Constant in equation 23 

Constant in equation 24- 

Constant in state equation (equation 13) 

Constant in state equation (equation 13) 

Discharge coefficient 

Specific heat at constant volume 

Ideal-gas specific heat at constant volume 

Critical- flow factor as defined in equation 1 6 

Perfect-gas critical-flow factor 

Constant in state equation (equation 13) 

Composition factor as defined in equation 22 

Mass flow rate per unit area 

Enthalpy (energy per unit mass) 

Molecular weight 

Mass flow rate 

Specific-heat ratio 

Number of components in mixture 

Pressure 

Specific gas constant 

Gao constant for calculating mass flow rate (equation 20) 

Gas constant for calculating volume flow rate (equation 21 ) 

Entropy 

Temperature 

Velocity 
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V Volume flow rate 

X Mole fraction of component in gas mixture 

Z Compressibility factor 

o: Constant in state equation (equation 13) 

Q^. Speed of sound in nozzle throat 

T Constant in state equation (equation 13 ) 

Constants in Ideal-gas specific-heat equation (equation 15 ) 
p Density 

Subscripts 

p Plenum 

t Nozzle throat 

i,j Species 

k Running index in equation (15) 

ANALYSIS 

In calculating the mass flow rate of a gas through a cri\ ical-fiow 
nozzle, certain assumptions are made. One is that the flow is isentropic 
from the plenum, where the gas is at rest, to the nozzle throat, where the 
flow velocity is sonic. Hie other assumption is that cne flow is one- 
dimensional. These assumptions are represented by the following equations: 


S p - S t “ 0 

u t 3 2 < H p - V 


( 1 ) 

( 2 ) 


Since the flow in the nozzle throat is sonic, an additional equation is: 



(3) 


Expressions for equations 1, 2, and 3 in terms of density and temperature 
will now be developed. 
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The differential expressions for entropy and enthalpy are given in 
reference 6 and are: 

dp (4) 

P 

dH = TdS + -4 dp + d® (5) 

p 2 w 

In terms of the compressibility factor, Z, equations 4 and 5 become 




where 


Z = Z(p,T) = 


Equations 6 and 7 are integrated along the path indicated in the following 


sketch. 
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In terms of these Integrations, equations 1 and 2 may be written as 



The term 3^ is the ideal-gas specific heat, and represents the specific 
heat the gas would have in the limit of zero density. 

Reference 7 gives an expression for the speed of sound in terms of 
the compressibility factor. Ibis can be used in equation 3 to yield: 
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In the calculation of the critical mass flow rate, the independent 
variables would be the plenum temperature and pressure. Equations 9 , 10, 
and 11 involve plenum temperature and density. For this reason, equation 8 
has to be solved implicitly for the plenum density. Once the plenum density 
is known, equations 9 > 10, and 11 can be solved for the nozzle throat 
density, temperature, and velocity. Hie iteration procedures used in 
solving these equations are the same as those used in reference 5* Once the 
gas properties in the nozzle throat are known, the mass flow rate per unit 
area can be calculated by the following equation. 


G t " ^ U t 


( 12 ) 


BASIC EQUATIONS 


One of the basic equations used in this calculation is the state equation 
that gives the compressibility factor as a function of density and temperature. 
For a one-component gas, there are many forms of this equation that can be 
chosen. Natural gas, however, is a mixture of many components. Furthermore, 
the composition of natural gas varies with time and region. It would 
therefore be convenient to be able to determine the coefficients of the state 
equation for a gas mixture from the coefficients of individual components. 

Such an equation is the one developed by Benedict, Webb, and Rubin in 
reference 8. This equation is as follows: 


Z - 1 + (B o - gj - % ) p + (b - ® T ) p2 + rt p5 + om 3 p ^ 1 


-rp 


(13) 


The coefficients in this equation are determined from the coefficients 
in the individual mixture components in the following manner (ref. 9). 
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Natural gas was assumed to consist of a mixture of seven components. 
These are methane, ethane, propane, 2-methyl propane, butane, nitrogen, 
and carbon dioxide. For pipeline gas, it would be unusual to find any 
appreciable amounts of hydrocarbons heavier than butane. The Benedict- 
Webb-Rubin coefficients for these gases are listed in Table I. 


Table I 

Benedict-Webb-Rubin Coefficients Fbr 
Components of Natural Gas 



Methane 

M*l6.043 

§Sg§ 

1 

Propane 

M=44.097 

2 -Methyl 

Propane 

M=58.124 

Butane 

M-58.124 

Nitrogen 

M=23.0l6 

Carbon 
Dioxide 
Mr 44. Oil 

-4 

Ax 10 
0 

.699525 

1.56707 

2.59154 

3.85874 

3.80296 

.449695 

1.03904 

3 o 

.682401 

1.00554 

1.55884 

2.20325 

1.99211 

.733661 

.790477 

C x 10 
0 

2.75763 

21.9427 

62.0993 

103.847 

121.305 

. 71953 ^ 

16.2657 

! a x 10"3 

2.98412 

20.8502 

57*248 

117.047 

113.705 

.900068 

3.94505 

b 

.867325 

2.85393 

5.77355 

10,889 

10.2636 

.508467 

.596615 

i—isa 

4.98106 

64.1314 

252.478 

559.777 

619.256 

1.07267 

24.6137 

a 

.511172 

1.00044 

2.49577 

4,41496 

4.52693 

1.19838 

3.23662 

r 

1.53961 

3.0279 

5.64525 

8,72447 

8.72447 

1.92451 



The above coefficients are for temperature in degrees Rankine and density 
in pound moles per cubic foot. Hie universal gas constant is taken as 10.7314 ft, ^ 
psi per °R per pound mole. The coefficients for all the gases except nitrogen 
and carbon dioxide are from reference 8. Those for nitrogen are from reference 
10 and those for cart* n dioxide were determined from a curve fit of the 
compressibility data in reference 3* 
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Another basic equation gives the ideal-gas specific heat as a function 
of temperature, flor a multicompound mixture, this specific heat is given by 


i R» 


where the ideal-gas specific heat of an individual component is given by 


+ 

R H o 


Za 


The values of are given in Table IX for the seven gases. 


Table II 

Constants in the Specific-Heat 
mial for Components of Natural Gas 
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These coefficients give the correct value of the ideal -gas specific heat 
over a temperature range of 360°R to 720°R. In all cases, these coefficients 
were determined from a curve fit of tabulated data, Hie methane data are 
from reference 12. These data represent a detailed analysis of spectroscopic 
data and agree with the data of the American Petroleum Institute (Reference 13 ) 
to within one percent. Hie ethane and propane data below 490°R are from 
reference 14. ethane and propane data above 490 °R are from reference 13 

as are the 2-methyl propane and butane data. The data for nitrogen are from 
reference 15 and the data for carbon dioxide are from reference 3* 

RESULTS AND DISCUSSION 

Hie results of these calculations are presented in terms of a dimensionless 
critical-flow factor which is defined as follows: 



If the gas were perfect, the expression for the critical- flow factor would 
only involve the specific-heat ratio and would be 


C* . 
perf 



(n+l)/(n-l) 


1/2 


(17) 


In terms of the critical-flow factor, the actual mass flow rate, m, of a 
gas through a critical- flow nozzle of geometric throat area, A^, is given by 


* - C D A t 



( 18 ) 


The upstream volume flow rate would be given by 


l 


12 


V P = C D A t 


C* Z 
P 



(19) 


In equations l8 and 19, the discharge coefficient mainly represents 
the effects of the non-isentropic and non-one-dimensional flow in the 
boundary layer of the nozzle. This coefficient is usually determined by a 
nozzle calibration and is generally plotted as a function of Reynolds 
number. Topical values of are between O.96 and 1.0. Hie results for 
the nozzle used in reference 16 indicate that is independent of Mach 

number from a value of 0.2, where the flow can be considered incompressible, 
to a value of 1, where the flow is critical and compressible. Since 
compressibility effects are negligible, real-gas effects should be negligible 
because, for the range of temperatures and pressures considered, the real -gas 
effects can be considered a correction to the compressibility effects. 

The specific gas constant in equation 18 is 


S 

m M 


( 20 ) 


where M is the molecular weight. This will give m in pounds per second 
if Pp is in pounds per sq. in., A^ is in square inches, and T^ is in 
degrees Rankine. 

The specific gas constant for equation 19 is 


r 

v M 


( 21 ) 


This will give V in cubic feet per second. Hie units of A. and T 
P t p 

are the same as those in equation 20. 

The critical-flow factor was calculated for methane and for nineteen 
natural-gas compositions. These compositions were determined from natural-gas 
samples taken from pipelines throughout the country. These data were supplied 
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in a private communication from Mr. Harry Schroyer of the American Meter 
Company. The range of composition of these gases is listed in Table III. 

Any alkane in the natural gas having more than 
four carbon atoms was ignored. The alkanes that were Ignored amounted to 
less than 0.4 percent by volume; in most cases, the amount was about 0.1 
percent by volume. 

Figure 1 presents the critical-flow factor for methane. Figures 2 and 
3 present the critical-flow factor for two typical natural-gas mixtures. 

The bulk of the natural gas found in pipelines has a composition similar 
to the composition of the natural gas of figure 2. The dotted portion of 
the 450°R curve on figure 3 represents a region in the calculation where 
the partial pressure of one of the components is above the saturation 
pressure. It is observed that in the limit of zero pressure, a 1^ percent 
variation of the critical- flow factor with temperature exists. This is due 
to the temperature variation of the ideal-gas specific heat. At a given 
temperature, the variation of the critical-flow factor with pressure is as 
much as 20 percent for the pressure range of 0 to 1000 pounds per square inch. 

Since the composition of natural gas is not fixed, and the critical-flow 
calculation is complex, it would be desirable to find a convenient empirical 
method to make this calculation. The method proposed in this paper makes use 
of a composition factor, f. For methane and the nineteen natural gas 
mixtures, both the critical-flow factor and the square root of the 
compressibility are linear functions of this composition factor. This factor 
is defined as follows: 

f = X CH + ~ +2X ch + 3X ch 

C 2 H 6 °°2 z C 3 H 8 C 4 H 10 

l 


( 22 ) 
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The X represents the mole fraction of the component whose chemical symbol 
appears as a subscript on the X. Hie accuracy of the critical-flow factor 
correlation was improved if the correlation was in terms of the product of 
the critical-flow factor and the square root of the compressibility factor 
rather than of the critical-flow factor alone. This correlation is 
represented by the following equation 

C* *Jz- a f + b (23) 

c c 

where a c and b Q are functions of plenum pressure and temperature and are 
listed in tables IV and V. In this equation, b c represents the contribution of 
the methane and a c f represents the contribution of the other components. (For 
a perfect gas with T = 4/3 , the value of b, would be O.6732). 

The coefficients in Eq. (22) and the form of the equation are entirely 
empirical. They result from an intuition regarding the relative contributions 
of the various components and several trial -and-error assumptions of numerical 
values of the coefficients. The values calculated by equation 23 agree with 
those calculated by the methods of this paper to better than 0.1 percent. 

In a similar manner, the square root of the compressibility factor can be 
represented by 


V Z = a f + b 

r 7 15 


(24) 


where a and b are functions of pressure and temperature and are listed 
z z 

in tables VI and VH. The values calculated by equation 24 agree with the values 
calculated by equation 13 to better than 0.5 percent. In most of the cases, 
the agreement is better than 0.1 percent. 

In terms of the factors calculated by equations 23 and 24, the 
expressions for mass flow rate and volume flow rate become 

(a„f + O P- 


m = 


Vt 


(ajf + 4,) 


and 


R T 
m p 


■ <Vt <\ f + W + V Vv p 


(25) 

( 26 ) 
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Since there appears to be no significant loss of accuracy in the use of 
equations 25 and 2b, their use is recommended rather than the more involved 
calculations described in the analysis section of this paper. 

In all the calculations in this paper, a state equation is involved. 

It should, therefore, be remembered that the accuracy of the results is 
probably no better than the accuracy of the state equation. 

SUMMARY 

. The mass flow rate of methane and nineteen natural gas mixtures through 
critical flow nozzles has been calculated. The calculation assumes the flow 
to be one-dimensional and isentropic. Hie pressure range is 0 to 1000 pounds 
per square inch and the temperature range is from 450 to 700 degrees Rankine. 

From a study of the results, a simple empirical method for making this 
mass flow rate calculation is proposed. This method would apply to natural 
gas mixtures whose composition is known and whose components have no more 
than four carbon atoms. 


TABLE III 

Range of Composition 
of Nineteen Natural Gases 


Gas 

Mole Fraction 

Methane 

Ethane 

Propane 

2 -Methyl Propane 
Butane 
Nitrogen 
Carbon Dioxide 

0.840 - O.967 
0.003 - 0.110 

Trace - 0.020 
Trace - 0.004 
Trace - 0.004 

0.001 - 0.023 
Trace - 0.017 
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-0.0456 

690 

-0.0)67 

-0.0)80 

-0.0)9) 

-9.0403 

-0.0416 

-0.0477 

-0.04)6 

-0.044) 

-0.0450 

-0.0434 

700 

-0.0167 

-0.0160 

-0.0)92 

-0.0404 

-0.0413 

-0.0423 

-0.04)4 

-0.0441 

-0.0447 

-0.0452 






• 

Ihhle V 










Value of b 

c 






Plenum 




Plenum praaaure. 

pela 





temperature, 











deg R 

0 

100 

200 

)0C 

400 

500 

600 

700 

800 

900 

4*0 

0.6719 

0.6715 

0.671) 

0.6712 

0.6713 

0.6717 

0.672) 

0.673) 

0.6747 

0.6767 

469 

0.6717 

0.6714 

0.6712 

0. 6712 

0.671) 

0.6717 

0.6724 

0.67)4 

0.6748 

0.6766 

4 70 

0.6714 

0.671? 

0.6711 

0 *711 

0.6714 

0.6718 

0.6723 

0.6734 

0.6747 

0.6764 

460 

0.6712 

0.6710 

0.6710 

0.67U 

0.671) 

0.6718 

0.6725 

0.67)4 

0.6747 

0.676) 

490 

0.4709 

0.6708 

0.6708 

0.6709 

0.671? 

0.6717 

0.6724 

0.6734 

0.6746 

0.6761 

500 

0.6707 

0.6706 

0.6706 

0.6708 

0.6711 

0.6716 

0.6723 

0.673) 

0.6743 

0.6759 

5 10 

0.6704 

0.6701 

0.6704 

0.6706 

0.6710 

0.6713 

0.6727 

0.6731 

0.674) 

0.6757 

* 70 

0.6701 

0.6701 

0.670? 

0.6704 

0.6708 

0.6714 

0.6771 

0.6730 

0.6741 

0.6735 

5)0 

0.6698 

0.6698 

0.6699 

0.670? 

0.6706 

0.6712 

0.6719 

0.6728 

0.67)9 

0.6732 

540 

0.6694 

0.6695 

0.6697 

0.6700 

0.6704 

0.6709 

0.6717 

0.6726 

0.67)6 

0.6749 

530 . 

0.6691 

0.6692 

0.6694 

0.6697 

0.6701 

0.6707 

0.6714 

0.6723 

0.67)4 

0.6746 

569 

0.6687 

0.6689 

0.6691 

0.6694 

0.6699 

0.6704 

0.6712 

0.6720 

0.67)1 

0.674) 

579 

0.6684 

0.6685 

0.6687 

0.6691 

0.6696 

0.6701 

0.6709 

0.6717 

0.6727 

0.67)9 

560 

0.6680 

0.6681 

0.6684 

0.6688 

0.6692 

0.6698 

0.6704 

0.6714 

0.6724 

0.67)3 

599 

0.667b 

0.6678 

0.6680 

0.6684 

0.6689 

0.6693 

0.6702 

0.6711 

0.6720 

0.67)2 

600 

0.6672 

0.6674 

0.6677 

0.6681 

0.6686 

0.6692 

0.6699 

0.6707 

0.6717 

0.6728 

610 

0.6668 

0.6670 

0.667) 

0.6677 

0.668? 

0.6688 

0.6693 

0.670) 

0.671) 

0.672) 

6?0 

0.666) 

0.6666 

0.6669 

0.667) 

0.6678 

0.6684 

0.6691 

0.6699 

0.6709 

0.6719 

6)0 

0.6659 

0.6662 

0.6665 

0.6669 

0.6674 

0.6680 

0.6687 

0.6693 

0.6705 

0.6715 

640 

0.6655 

0.6657 

0.6661 

0.6663 

0.6670 

0.6676 

0.668) 

0.6691 

0.6700 

0.6710 

650 

0.6630 

0.6651 

0.6636 

0.6661 

0.66G6 

0.6672 

0.6679 

0.6687 

0.6696 

0.6706 

660 

0.6646 

0.6649 

0.6652 

0.6657 

0.6662 

0.6668 

0.6673 

0.6683 

0.6691 

0.6701 

670 

0.6641 

0.6644 

0.6648 

0.6632 

0.6638 

0.6664 

0.6671 

0.6678 

0.6687 

0.6696 

680 

0.66)7 

0.6640 

0.664) 

0.6648 

0.665) 

0.6659 

0.6666 

0.6674 

0.668? 

0.6692 

690 

0.66)2 

3.66)5 

0.66)9 

0.6644 

0.6649 

0.6633 

0.666? 

0.6669 

0.6678 

0.6687 

700 

0.6677 

0.66)1 

0.66)3 

0.6619 

0.6644 

0.6651 

0.6657 

0.6663 

0.667) 

0.668? 


I 


1000 

-0.070b 

-0.0790 

-0.0)4) 

-0.0)76 

-0.0)99 

-0.0416 

-0.0477 

-0.04)6 

-0.044) 

-0.0446 

-0.0432 

-0.0433 

-0.0438 

-0.0460 

-0.0462 

-0.046) 

-0.0466 

-0.Q464 

-0.0464 

-0.0464 

-0.046) 

-0.0462 

-0.0461 

-0.0439 

-0.0437 

-0.0433 


1000 

0.6791 

0.6789 

0.6766 

0.676) 

0.6760 

0.6777 

0.6774 

0.6771 

0.6767 

0.6764 

0.6760 

0.6736 

0.673) 

0.6746 

0.6744 

0.6740 

0.67)3 

0.67)1 

0.6726 

0.6721 

0.6717 

0.6712 

0.6707 

0.6702 

0.6697 

0.6692 


Tibi® VI 


Valu* of ® z 


Plenum 
t- -j- rature, 


at* R 

0 

too 

100 

*50 

0 . 

-0.0151 

-0.09)0 

♦ 60 

0 . 

-0.011* 

-0.0*50 

♦ TO 

0 . 

-0.0110 

-0.0*5* 

♦ 00 

0 . 

-0.0101 

-0.0*11 

♦ ♦0 

0 . 

-0.0150 

-0.0)5) 

500 

0 . 

-0.0170 

-0.0)6* 

510 

0 . 

- 0 . 01 ** 

-0.01*1 

920 

0 . 

-0.015* 

-0.0)11 

5)0 

0 . 

-0.01*7 

-0.0)01 

5*0 

0 . 

-0.01)0 

-0.010) 

550 

0 . 

-0.01)0 

-0.01*6 

5*0 

0 . 

- 0 . 011 ) 

•0.0151 

5 TO 

0 . 

-0.011* 

-0.01)6 

500 

0 . 

-0.0110 

-0.011) 

550 

0 . 

-0.010* 

-0.0111 

*00 

0 . 

-0.0055 

-0.0100 

*10 

0 . 

-0.005* 

-0.0150 

*70 

0 . 

-0.0005 

-0.0100 

*)0 

0 . 

-0.0005 

-0.0171 

*♦0 

0 . 

-0.0001 

-0.01*1 

*50 

0 . 

-0.0077 

-0.019* 

**0 

0 . 

-0.007) 

-0.01*7 

*70 

0 . 

-0.0070 

-0.01*0 

*00 

0 . 

-0.00*7 

-0.01 )* 

*50 

0 . 

-0.006* 

-0.0110 

700 

0 . 

-0.00*1 

-0.0111 


Plenum 

temperature. 


de* R 

0 

100 

200 

*50 

1.0000 

0.5051 

0.5700 

*60 

1.0000 

0.5055 

0.575* 

*70 

1.0000 

0.550* 

0.5010 

*00 

l.CCOO 

0.5511 

0.502* 

♦ 59 

1.0000 

0.5510 

0.50)6 

500 

1.0000 

0.552* 

0.50*7 

510 

1.0000 

0.5525 

0.5*97 

510 

1.0000 

0.55)) 

0.50** 

5)0 

1.0000 

0.55)7 

0.5079 

5*0 

1.0C00 

0.55*1 

0.500) 

950 

1.0000 

0.55*5 

0.5051 

9*0 

1.0000 

0.55*5 

0.5050 

570 

1.0000 

0.5552 

0.550* 

900 

1.0000 

0.5555 

0.5510 

550 

1.0000 

0.5590 

0.551* 

*00 

1.0000 

0.55*0 

0.5521 

*10 

l.OCOO 

0.55*5 

0.552* 

*10 

1.0000 

0.55*5 

0.55)1 

*90 

1.0000 

0.55*7 

0.55)5 

6*0 

1.0000 

0.55*5 

0.55)5 

*90 

1.0000 

0.5571 

0.554) 

6*0 

1.0000 

0.557) 

0.55*7 

*70 

1.0000 

0.5579 

0.5590 

*00 

1.0000 

0.557* 

0.559) 

*50 

1.0000 

0.5570 

0.559* 

700 

1.0000 

0.5575 

0.5595 


Plenum preaaure, pale 


)00 

*00 

900 

*00 

-0.00)7 

-0.1175 

-0.19*1 

-0.1500 

-0.0770 

-0.1070 

•0.1*17 

-0.1750 

-0.0710 

-0.0505 

-0.125) 

-0.1*22 

•0.0*97 

-0.051 1 

•0.110* 

-0.1*70 

-0.0*10 

-0.00*2 

-0.1005 

-0.1)91 

-0.09*7 

-0.0700 

-0.1009 

-0.12*1 

-0.0520 

-0.072* 

•0.05)1 

-0.11*7 

-0.0*5) 

-0.0*79 

-0.000* 

-0.10*1 

-0.0*61 

-0.0*)0 

•0.000* 

-0.050* 

-0.0*)) 

-0.0905 

-0.0790 

-0.0916 

-0.0*0* 

-0.0991 

•0.0701 

-0.009* 

-0.0)02 

-0.0910 

•0.0*9* 

-0.079* 

-0.09*0 

-0.0*07 

•0.0*10 

-0.07*0 

-0.0))5 

-0.0*90 

-0.0975 

-0.0701 

-0.0)20 

-0.0*)2 

•0.09*9 

-0.0*99 

-0.0)0) 

-0.0*00 

•0.091* 

-0.0*21 

-0.0207 

-0.0)09 

•0.0*09 

-0.0909 

-0.0272 

-0.09*9 

•0.0*90 

-0.0552 

-0.0290 

-0.0)** 

-0.0*9* 

-0.0922 

-0.02*9 

-0.0)20 

•0.0*11 

-0.0*9* 

-0.02)) 

-0.0)11 

•0.0)50 

-0.0*00 

-0.0221 

-0.025* 

•0.0)70 

-0.0*** 

-0.0211 

-0.0201 

•0.0192 

-0.0*22 

-0.0201 

-0.02*0 

•0.0)9* 

-0.0*01 

-0.0152 

-0.0299 

•0.0)10 

-0.0)01 

-0.010) 

-0.02*9 

•0.0)0) 

-0.010) 



Tibi, vn 



Value of 



Plenum preaaure, 

pala 


)00 

*00 

900 

*00 

0.5*67 

0.9552 

0.9*)* 

0.9)15 

0.9*52 

0.9909 

0.9*70 

0.5)70 

0.571* 

0.901* 

0.9910 

0.9*19 

0.971* 

0.9*** 

0.995) 

0.9*0) 

0.9791 

0.90*9 

0.95*0 

0.990) 

0.5770 

0.909) 

0.9*10 

0.99K 

0.5705 

0.971* 

0.9**) 

0.9571 

0.5000 

0.971* 

0.9*00 

0.9*0) 

0.901) 

0.9752 

0.9091 

0.9*11 

0.502* 

0.97*0 

0.9712 

0.9*97 

0.90)7 

0.970* 

0.97)2 

0.9*01 

0.90*0 

0.979* 

0.9750 

0.9702 

0.9097 

0.9011 

0.9760 

0.972) 

0.90*7 

0.902* 

0.9702 

0.97*1 

0.5*75 

0.90)5 

0.9797 

0.9759 

0.500) 

0.90*0 

0.9*10 

0.9775 

0.9091 

0.909* 

0.9*2) 

0.9790 

0.9090 

0.90*9 

0.9*1* 

0.900* 

0.990* 

0.907* 

0.90*9 

0.9017 

0.9910 

0.9*02 

0.9090 

0.90)0 

0.991* 

0.9*90 

0.90*9 

0.90*1 

0.9922 

0.9*97 

0.907* 

0.9*92 

0.5927 

0.950* 

0.900) 

0.90*2 

0.99)1 

0.9911 

0.9091 

0.9072 

0.99)0 

0.9917 

0.9090 

0.9001 

0.99*0 

0.9922 

0.9909 

0.9009 


700 

*00 

900 

1000 

-0.2*0* 

-0.2991 

-0. )5** 

• 0.* 1*2 

-0.2199 

-0.2*** 

-0.1119 

-0.1610 

-0.197* 

-0.2)60 

•0.27*2 

-0.1175 

-0.1791 

-0.212) 

-0.2*69 

-0.2021 

-0.10)1 

-0.1921 

-0.2225 

-0.2512 

-0.1*9) 

-0.1752 

•0.2010 

-0.2200 

-0.1)71 

-0.1*0* 

-0.10*1 

'0.2079 

-0.1209 

-0.1*7* 

-0.1*07 

-0.1900 

-0.1170 

-0.1 )*0 

-0.1552 

-0. 17*3 

-0.100* 

-0.1259 

-0. 1*1) 

-0.1*06 

-0.1010 

-0.11*9 

-0.1)27 

-0.1*0* 

-0.09*2 

-0.1000 

-0.12)1 

-0.1)77 

-0.0001 

-0.1015 

-0.11*9 

-0.1201 

-0.0029 

-0.09*9 

-0.1072 

•0.119* 

-0.077* 

-0.0009 

-0.100) 

•0.111* 

-0.0720 

-0.0*15 

-0.09*1 

-0.10*5 

-0.0*09 

-0.07*5 

-0.0001 

-0.0900 

-0.0*** 

-0.07)9 

-0.00)1 

-0.0921 

-0.0*10 

-0.0097 

-0.97*1 

-0.00*7 

-0.0577 

-0.0*50 

-0 17)9 

-0.081 B 

-0.09** 

-0.0*2) 

-0.0*90 

-0.0772 

-0.0917 

-0.0909 

-0.0*60 

-0.07)0 

-0.0*91 

-0.0959 

-0.0625 

-0.0*91 

-0.0*** 

-0.09)0 

-0.0591 

-0.065* 

-0.0**) 

-0.090) 

-0.096) 

-0.0*21 

-0.0*21 

-0.0*79 

-0.0515 

-0.05B9 


700 

000 

900 

1000 

0.9195 

0.9075 

0.0959 

0.00)7 

0.9261 

0.9192 

0.90** 

0.09)0 

0.9)20 

0.9221 

0.912* 

0.9020 

0.917) 

0.92*1 

0.9195 

0.9109 

0.9*20 

0.91)9 

0.9259 

0.9101 

0.9*6* 

0.9109 

0.9)1* 

0.92*5 

0.990) 

0.9*15 

0.91*9 

0.910* 

0.9919 

0.9*77 

0.9*1* 

0.9)57 

0.9572 

0.9915 

0.9*59 

0.9*0* 

0.9*0) 

0.9590 

0.9*99 

0.9*50 

0.9*11 

0.9902 

0.95)9 

0.9*90 

0.9*5* 

0.9*12 

0.95*9 

0.9527 

0.9*00 

0.96)9 

0.9599 

0.99*2 

0.9702 

0.9*6* 

0.9*20 

0.959) 

0.972) 

0.9*00 

0.905* 

0.9*22 

0.97*1 

0.9709 

0.9070 

0.9*49 

0.9799 

0.9729 

0.9701 

0.9*7* 

0.9776 

0.97*0 

0.9722 

0.9*90 

0.9791 

0.9766 

0.97*2 

0.9719 

0.9009 

0.9702 

0.9760 

0.97*0 

0.9019 

0.9797 

0.9777 

0.9750 

0.90)1 

0.9012 

0.9791 

0.977* 

0.90*1 

0.9025 

0.9000 

0.979) 

0.909* 

0.90)0 

0.9022 

0.900* 

0.9009 

0.90*9 

0.9019 

0.902) 

0.907* 

0.90*1 

0.90*0 

0.90)6 


/ 
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PLENUM PRESSURE, PSI t 
Figure 1. - Critical flow factor for methane. 


PLENUM PRESSURE, PSIA 

Figure 2. - Critical flow factor for a natural gas. 
Composition by volume is: 0.9535 methane. 
0.0296 ethane, 0.0046 propane, 0.0007 2-methyl 
propane, 0.0006 butane, 0.0040 nitrogen, and 
0. 0070 carbon dioxide. 


, PLENUM 
TEMPEF- 
! ATURE 
| °R 
i 450 


CRITICAL .76 
ROW 
FACTOR, 

C* • 


PLENUM PRESSURE, PSIA 

Figure 3. - Critical flow factor for a natural gas. 
Composition by volume is.- 0. 8850 methane, 
0.0795 ethane, 0.0110 propane, 0.0007 2-methyl 
propane, 0.0017 butane, 0.0221 nitrogen, and 
0 carbon dioxide. 







